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Abstract. The gross features in the distribution of stars as well as warm (T S 100 K) interstellar dust in the Galactic disk have 
been investigated using the recent mid infrared survey by Midcourse Space Experiment (MSX) at 8, 12, 14 & 21 pm bands. 

An attempt has been made to determine the location of the Galactic mid-plane at various longitudes, using two approaches : (i) 
fitting exponential functions to the latitude profiles and (ii) statistical indicators. The former method is successful for the inner 
Galaxy (-90° < I < 90°), and quantifies characteristic angular scales, y, along latitude. These ys have been translated to linear 
scale heights (zh) and radial length scales (R;) using geometric description of the Galactic disk. The distribution of warm dust 
in the Galactic disk is found to be characterised by R/ < 6 kpc and 60 100 pc, consistent with other studies. 

The location of the Galactic mid-plane as a function of longitude (in all 4 MSX bands), for stars as well as warm dust, has been 
searched for signatures of warp-like feature in their distribution, by fitting sinusoid with phase and amplitude as parameters. 

In every case, the warp signature has been detected. Carrying out an identical analysis of the DIRBE/COBE data (with lower 
angular resolution) in all its ten bands covering the entire infrared spectrum (1.25-240 pm), also leads to detection of warp 
signatures with very similar phase as found from the MSX data. Our results have been compared with those from other studies. 
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1. Introduction 

The study of structural features of the Galactic disk is useful 
in understanding the spatial distribution of its constituents, viz. 
stars, interstellar gas and interstellar dust. Our Galaxy consists 
of the following kinematically distinct components: spheroidal 
(bulge, halo) and disk components (thin disk, thick disk). The 
disk component of the Galaxy has been studied extensively us¬ 
ing several infrared surveys. Kent et al. lfT99Tt obtained the 
three dimensional distribution of our Galaxy based on data at 
2.4 pm from the Spacelab IRT. Freudenreich J1996t used the 
Diffuse Infrared Background Experiment (DIRBE/COBE) data 
in J, K, L and M bands to model the Galactic disk. 

One of the subtler features of the Galaxy is that of the 
warp structure, first detected through the 21 cm line emission 
of H I (Burke [11)571 KeiT 11957b . Various theories have been 
put forward to explain this phenomenon : gravitational interac¬ 
tion with a companion galaxy (Burke n~957l Kerr i 19571 Hunter 
& Toomre ri969l and others), effect of massive halo (Binney 
[I9781Q991, misaligned halo-disc axes leading to cosmic in¬ 
fall (Ostriker & Binney 179891 Jiang & Binnev IT999t . inter- 
galactic magnetic fields (Battaner et al. rr^rrwTb . dynamical 
friction between misaligned rotating halo and disc (Debattista 
& Sellwood ll999> . and more recently, intergalactic accretion 
flows onto the disc (Lopez-Corredoira et al. I2002at . Porcel 
et al. <1997t have explained the Galactic stellar warp using 
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the tilted ring model of Rogstad et al. ( fD7?b . The warp in 
the stellar disk was detected from the distribution of sources 
in the Infrared Astronomical Satellite (IRAS) Point Source 
Catalog (Djorgovski & Sosin ll9^ . Using the DIRBE data, 
Ereudenreich et al. <1994> . found evidence of warp for the stel¬ 
lar and the interstellar dust components of the Galaxy. The 
old stellar populations near the Galactic plane have been anal¬ 
ysed for Galactic flare and warp from the Two Micron All Sky 
Survey (2MASS) by Lopez-Corredoira et al. J2002bll . 

Recently, the Midcourse Space Experiment (MSX) sur¬ 
veyed the entire Galactic plane in four mid infrared spectral 
bands (8-21 pm) with high (~18") angular resolution (Price et 
al. lMTTb . Taking advantage of the complete longitude coverage 
and superior angular resolution of this new dataset, we have at¬ 
tempted to extract information about structure of the Galactic 
disk and search for possible warp signature from the distribu¬ 
tion of stellar as well as warm interstellar dust components of 
our Galaxy. Eor a direct and meaningful comparison, an iden¬ 
tical analysis has also been carried out with the DIRBE data, 
albeit with poorer angular resolution. 

In Section 2, we discuss the MSX and DIRBE datasets used 
for the present study. The data analysis procedures followed to 
extract structural features of the Galaxy and possible signatures 
of warp, have been described along with the results in Section 
3. Based on a geometric description of the Galactic disk, a rela¬ 
tion between the linear and the angular scale heights have been 
established in Section 4. Our results have been compared with 
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others’ in Section 5. In Section 6, a brief summary has been 
presented. 

2. Data Sets 

The main objectives of the present study are to quantify param¬ 
eters that characterise the structure of Galactic disk and search 
for possible signatures of warp from the distributions of stellar 
and warm interstellar dust components. The MSX data prod¬ 
ucts have been used extensively. In addition, DIRBE data have 
also been used for comparison. 

2.1. Midcourse Space Experiment (MSX) 

The MSX carried out mid infrared observations of the Galactic 
plane (|b| < 5°) in four bands A, C, D and E with /l(A/l) corre¬ 
sponding to 8.28 (3.36), 12.13 (1.72), 14.65 (2.23), and 21.34 
(6.24) jum, respectively (Price et al. I2001t . While the MSX 
Point Source Catalog, MSX PSC (Egan et al. I1999t . has been 
used to study the distribution of point sources, the panoramic 
images in the four bands (pixel size 36"x36") have been used 
to study the diffuse emission. These images have already been 
corrected for the zodiacal emission using an updated version 
of CBZODY model described in Kennealy et al. J1993t . The 
MSX PSC and the panoramic images were taken from IPAC 
(http://irsa.ipac.caltech.edu/applications/MSX/l 
and the MSX team. 

We selected 4 samples of sources from the MSX PSC (one 
for each band), each based on the following criteria : (1) posi¬ 
tive detection of the source in the particular band (quality flag 
> 2); (2) source located within the latitude range |b| < 2.75°, 
in order to avoid regions with non-uniform sky coverages; and 
(3) sources brighter than the completeness limit of the cata¬ 
logue for that band, viz. 0.158 Jy, 1.58 Jy, 1.58 Jy and 3.16 
Jy for 8, 12, 14 and 21 pm bands, respectively (as determined 
from differential logN-logS distributions). The completeness 
limit in A band corresponds to 6.42 magnitude at 8 pm, us¬ 
ing the magnitude scale given by Egan et al. (I1999t . Our sam¬ 
ples consist of 164723, 25467, 18239 and 14254 sources in 8, 
12, 14 and 21 pm bands, respectively. The 8 pm band of MSX 
PSC being the most sensitive, the corresponding sample is the 
largest. Using MSX PSC data, Lumsden et al. Il2()()2> have plot¬ 
ted known stars, planetary nebulae, H II regions and young 
stellar objects (YSOs) on colour-colour diagrams in Eigures 

3, 4 & 5 of their paper. We have used a colour cut of the 
form logio(T’i4/T'i2) = 0.87 logio(T' 2 i/T’ 8 ) - 0.61 (equation of 
a straight line on a log-log plot; stars are located to the left of 
this line) which separates the HII regions, etc from the stars, as 
evident from their Eigure 3 (right panel). Using this colour cut 
we find the fraction of non-stellar objects to be ~ 13%. This 
is similar to IRAS Point Source Catalog, where ~ 15% of the 
sources detected in both the mid infrared bands (12 & 25 pm) 
are non-stellar objects (Beichman et al. IMH}. 

2.2. Emission from DIRBE Sky Survey 

The DIRBE instrument onboard COBE spacecraft carried 
out a full sky survey in 10 photometric bands centered at 


1.25, 2.2, 3.5, 4.9, 12, 25, 60, 100, 140 and 240 pm, albeit 
with limited angular resolution of ~ 0.7° (Boggess et al. 

[7^ . We have used the Zodi-subtracted Mission Average 
(ZSMA) maps here, which were obtained from LAMBDA 
(bttp: //lambda. gsf c. nasa. gov/product/cobe/dirbejtroducts. cfn 

3. Analysis and Results 

3.1. MSX Point Source Catalog (PSC) Sources 

3.1.1. Estimated distribution of stellar types & distance 
for our sample 

In order to estimate the distribution of stellar types of the mid 
infrared sources in our sample and their distances, we have 
used a modified version of the Besan 5 on model of stellar popu¬ 
lation synthesis (A. Robin, private communication). An earlier 
version of the model has been described in Robin et al. (| 2003 . 

The modified version of the model includes predictions of star 
counts in two mid infrared bands corresponding to the LW2 
(7 pm) & LW3 (15 pm) filters of ISOCAM, which are crucial 
for the present study. This model synthesizes stars belonging 
to the populations of the thin disk, the thick disk, the stellar 
halo and the outer bulge including the warp and flare features. 

This model has also been used to estimate the contribution of 
all stars, including those not detected by MSX PSC, to the total 
mid IR emission (see section 3.2.1). Simulations have been car¬ 
ried out over 10x10 deg^ area for |b| < 5° and 30° < / < 330° 
on a grid with 30° spacing. The simulations covered distances 
upto 20 kpc from the Sun and stars upto apparent magnitude of 
14.0 at 7 pm. The extinction due to interstellar dust has been 
neglected in the model simulations. 

As mentioned earlier, the completeness limit of MSX PSC 
for A band (8 pm) corresponds to 6.4 magnitude at 8pm, [8]. 

This translates to a value of 6.5 magnitude at 7pm, [7], for the 
stellar types in our sample (the mean difference of the ([7] - 
[8]) colour, is ~ 0.1 mag; for various effective temperatures 
covered in our sample, as evident below). Accordingly, we have 
considered the predicted star counts from the simulations with 
[7] < 6.5, to be representative of our MSX PSC sample. 

The mean Galactocentric distance for our sample, is found 
to vary from ~ 4 kpc towards the Galactic centre (/ ~ 0°) to ~ 

11 kpc towards the anti-centre direction (I ~ 180°). Most of the 
sources are dominated by K and M giants (31% K & 59% M 
giants). The rest are AGB stars (6%) and earlier than G dwarfs 
(4%). 


3.1.2. Latitude distribution of stellar sources and warp 
signature 

The structure of the Galactic disk has been studied using the 
local distribution of stellar sources of our sample by construct¬ 
ing latitude profiles (LP) of stars for 36 longitude bins, each 
10° wide. The bin size along latitude is 0.1°. The stellar den¬ 
sity distribution in the Galaxy is generally represented by an 
exponential function of the form n(z, R) oc e •* , where z 

is the height above the Galactic plane, R is the Galactocentric 
distance, Zh and Ri represent the characteristic scale height and 
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radial scale length, respectively. However, all relevant measure¬ 
ments are (approximately) heliocentric, and conveniently pre¬ 
sented in Galactic coordinates (/, b). Later (Section 4), we shall 
relate the linear and the angular scales using a geometric de¬ 
scription of the Galactic disk and show that exponential func¬ 
tions provide good hts to the latitude prohles too. Let N(b) be 
the observed latitude prohle representing the number of stars as 
a function of latitude, b, for a given longitude zone. At hrst, we 
explore both the exponential and gaussian functions for htting 
N(b), as : 

N(b) - UeC ~ +6e + eeb (1) 

Nib) - UgC ^ (2) 

where a, (3, y, 5 and e are the free parameters. The subscripts 
e and g refer to the exponential and gaussian functions, respec¬ 
tively. The parameters Ue & ag represent the amplitude; j3e & 
j3g represent the position of the Galactic midplane dehned by 
the MSX PSC stars; je & jg represent the angular scale height; 
6e, 5g and e^,, eg represent the background, if any, limited to 
linear dependence on latitude. Generally, /3 - 0 is expected if 
the peak of the latitude prohle appears at h = 0. Any non-zero 
value signihes a shift. The httings have been carried out for all 
the 36 longitude bins using a non-linear least squares method. 
For each of the four wavebands, good hts were obtained for 
the inner half of the Galactic plane towards the Galactic center 
i.e. -90° < / < 90° and very poor hts for the outer half, i.e. 
towards the Galactic anti-center 90° < I < 270°. Predictions 
of the modihed Besan^on model also show the same trend for 
the distribution of stars ([7] < 6.5 mag). Figure ^ shows the 
latitude prohles (MSX PSC, A band) for three representative 
longitude bins (30°, 60° & 180°), including functional hts for 
the two longitudes in the inner Galaxy. For all the four bands, 
the exponential function hts the latitude prohles better than the 
gaussian function, which is consistent with our understanding 
of the Galactic disk. Consequently, the results pertaining to the 
exponential function only have been presented here. Table ^ 
lists the mean value and standard deviation of the parameter 
y for each waveband, as obtained from the 18 longitude bins 
covering the inner Galactic plane. 

Next we explore possible signatures of warping in the 
Galactic disk. To quantify any evidence for the Galactic warp 
from the local distribution of stars, we search for systematic 
variation of fis with longitude. Djorgovski & Sosin ( I1989t . used 
a sinusoid function to ht the distribution of IRAS point sources. 
A similar approach has been adopted here, but without the ad¬ 
ditive constant term, thereby assuming that the Sun is located at 
the local midplane. Thus, the function used to ht the longitude 
variation of the prohle peak is : 

j8(0 = ySo X sinil - 0) (3) 

where fi represents the best ht to Peil) of eqn (1) obtained for 
-90° < / < 90°. Although useful, this is admittedly an approx¬ 
imate approach due to the fact that the measurements are helio¬ 
centric (hence, lines of sight along different longitudes cover 
different galactocentric distances of the disc), whereas a true 


Table 1. Characteristic angular scale, y, of latitude prohles ob¬ 
tained by htting exponentials to the MSX PSC sample, MSX 
emission and DIRBE emission for the 18 longitude bins cover¬ 
ing -90° < / < 90°. 


Wavelength 


Exponential 

jim (Band) 


yg in deg(s.d.) 


MSX PSC 

MSX Diffuse emission 

8.3 (A) 

1.29(0.67) 

1.69(0.45) 

12.1 (C) 

1.27(0.57) 

0.92(0.41) 

14.7 (D) 

1.25(0.70) 

0.71(0.22) 

21.3 (E) 

1.04(0.96) 

0.56(0.28) 

DIRBE Emission 

1.25 


1.94(1.18) 

2.2 


2.74(0.78) 

3.5 


1.93(0.54) 

4.9 


1.57(0.47) 

12 


1.29(0.67) 

25 


1.00(0.46) 

60 


0.92(0.33) 

100 


1.23(0.77) 

140 


1.19(0.69) 

240 


1.30(0.87) 


sinusoid due to any warp feature is expected only for a galac¬ 
tocentric observer. We crudely relate the parameters Po and (ft 
to quantify the local warp signature (hereafter WS) indicating 
amplitude and phase respectively. This phase, cp, corresponds 
to the position angle (longitude) of the line of intersection be¬ 
tween the b-Q plane and the warped Galactic plane. In order 
to carry out a quantitative comparison with a ‘null’ hypothe¬ 
sis (i.e. WS absent), we also ht a constant function, pH) - k, 
to the same data sets. The resulting best ht parameters along 
with their standard errors have been presented in Table ElThe 
data and the htted functions are displayed in Figure E The WS 
has been clearly detected in all the 4 bands as is evident from 
the ratio of;^f^s per degree of freedom {R > 1). On repeating 
our analysis with a larger sample of MSX PSC sources with 
latitude cut relaxed to \b\ < 5°, again we detect WS with very 
similar parameters. 

Since only the inner half of the Galactic plane could be 
studied using the above approach, we pursue additional mea¬ 
sures of statistical nature to quantify the location of the Galactic 
midplane, which could be effective for the entire longitude 
range. We have explored two statistical measures of central ten¬ 
dency, viz., the mean and the median latitude. The mean refers 
to the average latitude of the MSX PSC sources of our sam¬ 
ple for a particular longitude bin. The median latitude refers 
to a value such that the total number of stars located at higher 
and lower latitudes are equal. The mean is a good measure of 
central tendency for roughly symmetric distributions but can 
be misleading for skewed distributions. The median should be 
more informative for such distributions. 

The mean and median latitudes for different Galactic lon¬ 
gitude bins have been displayed as histograms in Figure 0 A 
systematic deviation from the nominal Galactic plane {b = 0) is 
clearly visible in this hgure. Both the statistical indicators (SI) 
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are positive over one half of the longitude range ~ 0° -190° and 
negative in the other half. The positive and the negative peaks 
occur at I ~ 90° - 100° and ~ 260° - 280°, respectively. This 
is the general trend for all the 4 bands of MSX. Similar to the 
case of p above for the exponential fits (eqn. 3), we fit a sinu¬ 
soid (with bo and (p as parameters) to the mean and the median 
latitudes giving equal weightage to each longitude bin. Once 
again, the ‘null’ hypothesis is also explored for comparison. 
The resulting best fit parameters along with standard errors, for 
all the cases have been presented in Table 0 The correspond¬ 
ing fitted functions for respective bands have been overplotted 
in Figure 0 The signature of warp is clearly indicated from 
both the statistical indicators (Sis) for all the 4 MSX bands. 
This is also evident from the ratio of the (R > 1). These 
findings are statistically further strengthened by the fact that 
the values of the phase, 0, are consistent (within errors) among 
all four bands of MSX and also matches that expected from 
the known features of the Galactic stellar warp (Djorgovski & 
Sosin [T98^ see Section 5.2.1). The SI approach seems to be 
more successful than LP. This could be because the statistical 
determinations of central tendencies are inherently more im¬ 
mune to fluctuations than procedures involving fitting of a spe¬ 
cific function. 

The amplitude of the sinusoid is lower for the SI than the 
jSs obtained from the fitting of the latitude profile. Among the 
mean and median, the displacement is larger for the latter case 
than the former. The mean provides a better quality fit to sinu¬ 
soid than the median, based on comparison of corresponding 
reduced ;^f^s. 

3.2. Diffuse emission from MSX images 

Diffuse mid infrared emission is expected from the warmer (T 
^ 100 K) interstellar dust heated by nearby luminous stars. 
The emission in the mid infrared (obtained from the MSX 
panoramic images) contains contributions from stellar sources 
also. In the next subsection, we estimate the contribution of 
point sources, which turns out to be a negligible fraction of the 
total emission detected. 

3.2.1. Contribution of point sources to the MSX data 

The panoramic images from MSX includes the contribution 
from the point sources, some of which have been detected and 
are catalogued in the MSX PSC. The sources from our MSX 
PSC sample contribute nearly 4.5%, 1.9%, 3.5% and 2.9% of 
the total mid infrared emission in the 8, 12, 14 and 21 /um 
bands, respectively. (The corresponding fractions considering 
all the MSX PSC sources are : 4.7%, 1.9%, 3.8% & 3.1%.) 
However, these correspond only to the sources brighter than the 
completeness limit of the catalog. In other words, there are (1) 
unresolved sources (MSX angular resolution is ~18") and (2) 
sources fainter than the sensitivity limit, contributing to the in¬ 
tensity in these images, apart from the cataloged point sources. 
To estimate the contribution of unresolved and faint stellar 
sources in the mid infrared MSX data, we have used the mod¬ 
ified Besangon model described earlier (Section 3.1.1). The 


contributions from stars in the selected latitude range for vari¬ 
ous longitude zones have been determined from the model sim¬ 
ulations. Integrated fluxes at 7 //m from all stars within a helio¬ 
centric distance of 20 kpc and apparent brightness higher than 
m-urn ([7] < m/im, where m;™ = 6.5, 7.0,...., 14.0), have been de¬ 
termined. It was noted that beyond m/,m=12, the increase in the 
value of the integrated flux was insignificant (~ 2%). The ra¬ 
tio of integrated fluxes for sources with [7]< 14.0 (representing 
the entire Galaxy) and [7]< 6.5 (MSX PSC completeness at A 
band; section 3.1.1) is found to be ~ 1.4. Carrying out similar 
model simulations at 15 /rm and considering the completeness 
of MSX PSC at D band, we And the ratio of corresponding in¬ 
tegrated fluxes is ~ 1.3. It may be noted that the values of these 
ratios are in fact upper limits, since the interstellar extinction 
has been neglected in the model simulations. The above implies 
that, the contribution of all stars towards the total emission is 
~ 6.3% & ~ 4.6% for 8 & 14 yum MSX bands respectively. 
If we adopt a value of 1.4 for the ratio of integrated fluxes for 
the 12 & 21 fim bands also, then the stellar contributions in the 
respective bands would be ~ 2.7% & ~ 4.1%. Therefore, we 
conclude that most of the total mid infrared emission detected 
in the MSX bands are of diffuse origin. From a study of the dif¬ 
fuse infrared emission (for |b| > 10°) of our Galaxy using the 
IRAS data, Boulanger & Perault (I1988t found that the point 
sources account for 10% + 2% of the 12 fim flux. 


3.2.2. Latitude distribution of diffuse emission and 
warp signature from MSX data 

The latitude distribution of diffuse emission in the 4 MSX 
bands, has been studied following an analysis similar to that 
used for the stars (Section 3.1.2), but for |b| < 5°. Exponential 
and gaussian functions (see Eqns (1) and (2)) have been fitted 
to the latitude profiles of intensity at respective bands. The fits 
are found to be good for the inner Galaxy (-90° < / < 90°) but 
very poor towards the outer Galaxy, similar to the stellar case. 
The latitude profiles of intensity in A band for three represen¬ 
tative longitude bins (10°, 60° & 180°) have been presented 
in Eigure 0] which also displays the fitted functions for the 
two bins in the inner Galaxy. Again, we And that the exponen¬ 
tial provides better fits than the gaussian for all the cases, as 
expected. Henceforth, we discuss the exponential only. The av¬ 
erage values (& standard deviation) of the angular scale height 
parameter, je, from the 18 longitude bins, have been listed in 
Table ^ It may be commented here that, the inferred je should 
be immune to any residuals from the zodiacal subtraction, be¬ 
cause : (i) the angular scale of zody emission is much larger 
than typical y/, & (ii) the background terms (constant & linear 
in b) should absorb all components varying slowly with b. 

In order to search for WS in the diffuse emission, we fitted 
sinusoids to the locations of profile peak at various longitudes, 
for all 4 MSX bands (similar to the stellar case; Eqn. 
3). Once again, the WS has been detected. The best fit values of 
the parameters of the WS, viz., amplitude and phase are listed 
in Table |2l The ySf'^^(/)s and corresponding best fitting sinu¬ 
soids are presented in Eigure |3 
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We have also pursued the statistical approach (SI; similar 
to the stellar case) to extend the search for WS over the entire 
Galactic longitude range. Two measures of central tendencies 
have been defined for the latitude distribution of intensity at 
various 1. The intensity weighted central latitude of the 

diffuse emission is defined by : 

5 ° 5 ° 

BintUD = X Z 

bi=-5° 

where bi and /j are the latitude and the diffuse emission inten¬ 
sity for the pixel i, considering all pixels within the particular 
longitude bin /. The median latitude, B,„eii{l), is defined to be 
that latitude, on both sides of which the areas under the inten¬ 
sity vs latitude curve are equal, 

B,„eM 5 " 

fiAb = ^ fjAb (5) 

where Ab is the pixel size along latitude (36")- This equation 
(5) is solved to determine Bmed for each of the 36 longitude 
bins. The statistical indicators (SI) for the Galactic midplane, 
Bintwt & Bmed, as a function of Galactic longitude for each MSX 
band are presented in Figure These Sis for diffuse emission 
have also been fitted with sinusoids to search for WS. 

Table 0 lists the best fit parameters corresponding to the 
extracted WS for the SI. We observe that even in this case of 
diffuse emission from warm interstellar dust, there is a clear 
warp-like signature in all the four bands. Interestingly, the val¬ 
ues of the phases are very similar to the stellar case. For both 
the Sis we find that C and E bands show similar amplitudes 
within errors. However, for the D band we find that the Bmej 
shows a larger amplitude of WS than and vice-versa for 
the A band. This is in contrast to the stellar case where the me¬ 
dian showed a larger amplitude of WS in all the 4 bands. For 
the case of SI of central tendencies, a better quality of fit to si¬ 
nusoid is obtained for Bmed than Bi„twt, for almost all the cases. 

3.3. Latitude distribution of emission and warp 
signature from DIRBE data 

We have searched for WS from the emission of the Galactic 
plane over the entire infrared waveband using the DIRBE im¬ 
ages following an analysis identical to that for MSX images 
(Section 3.2.2). Despite poorer angular resolution, the large 
number of DIRBE wavebands (10) covering near (1.25 pm) 
to far (240 pm) infrared makes it a good candidate for compar¬ 
ison. The extracted characteristic angular scales obtained from 
fitting latitude profiles (LP) with exponential functions, y^,, for 
the inner Galaxy have been listed in Table ^ The Sis for the 
Galactic midplane, Bi„t„, & B^ed for the entire longitude range, 
for each DIRBE band have been presented in Figure 0 These 
Sis have also been searched for WS and the resulting best fit¬ 
ting sinusoids have been overplotted on the same figure. 

The WS has been detected in all the three (one LP and two 
SI) cases in each of the DIRBE bands. Even the phases of the 
WS, as evident from Table |2l(LP) (& 0 SI), obtained for 9 (8) 
of the 10 bands are consistent with those obtained from MSX 


Table 2. The values of the best fit parameters of the warp sig¬ 
nature, WS (sinusoid with amplitude /5o & phase <p) to the lon¬ 
gitude variation (-90° < I < 90°) of locations of the peaks of 
latitude profile, LP, obtained from the exponential fitting to : 
the MSX PSC sources (|b| < 2.75°), diffuse emission in MSX 
bands (|b| < 5°), and emission in DIRBE bands (|b| < 5°). The 
free parameter, k, for the null hypothesis is also tabulated. The 
ratio of per degree of freedom, R, for the best fitting null 
case and the WS case, is also presented. A value of R> 1 indi¬ 
cates that the sinusoid is a better fit. 


A 

Bo 

<!> 

k 

R 

(pm) 

(deg) 

(deg) 

(deg) 


MSX PSC 

8.3 

0.44 ±0.11 

3.1 ± 14 

-0.04 ±0.10 

1.95 

12.1 

0.39 ±0.10 

3.5 ± 16 

-0.00 ±0.10 

1.80 

14.7 

0.47 ±0.15 

3.2 ± 19 

-0.04 ±0.13 

1.53 

21.3 

0.42 ±0.13 

6.7 ± 18 

-0.10±0.11 

1.53 

MSX Diffuse Emission 

8.3 

0.40 ±0.10 

18 ± 14 

-0.13 ±0.09 

1.81 

12.1 

0.57 ±0.16 

11 ± 17 

-0.15 ±0.14 

1.58 

14.7 

0.44 ±0.14 

6.2 ± 18 

-0.13 ±0.12 

1.60 

21.3 

0.69 ±0.19 

-6.3 ± 16 

0.07 ±0.17 

1.75 

DIRBE Emission 

1.25 

1.33 ±0.37 

72± 16 

-0.71 ±0.29 

1.30 

2.2 

0.51 ±0.08 

19 ±9.5 

-0.08 ±0.10 

3.11 

3.5 

0.38 ± 0.06 

15 ±9.4 

-0.04 ± 0.07 

3.20 

4.9 

0.34 ± 0.07 

15 ± 11 

-0.04 ± 0.07 

2.63 

12 

0.36 ±0.11 

18 ± 17 

0.12 ±0.09 

1.50 

25 

0.39 ±0.16 

11 ±23 

-0.06 ±0.12 

1.30 

60 

0.45 ±0.15 

12± 19 

-0.07 ±0.12 

1.46 

100 

0.40 ±0.11 

13± 16 

-0.10 ±0.10 

1.65 

140 

0.37 ±0.11 

16± 17 

-0.12 ±0.09 

1.52 

240 

0.37 ±0.11 

16± 17 

-0.12 ±0.09 

1.48 


data for the stars as well as the diffuse mid infrared emission. 
It may be noted that at 1.25 pm, the only band where the WS 
phase (& also the amplitude) deviates significantly from the 
other bands, the errors are the largest. This could be because 
the emission at 1.25 pm originates mainly from the stars, which 
would be affected most by extinction due to interstellar dust. 

4. Characteristic anguiar and iinear scaie heights 

In the present study, we have determined the characteristic an¬ 
gular scales of emission along the Galactic latitude, which need 
to be translated to linear scales for drawing physical inferences. 
Eor this purpose, we have carried out calculations based on ge¬ 
ometric description of the Galactic disk. The Galactic disk has 
been considered to be of radius 14 kpc and thickness 700 pc, 
with Rq = 8 kpc. Our description neglects the effects of the 
Galactic bulge, the halo as well as any warp feature. The effect 
of these simplifications on our results is expected to be negligi¬ 
ble. 

We need to estimate the integrated emission as a function 
of Galactic latitude for a given longitude, I{b), with linear scale 
height (zh) and radial scale length (Ri) as input parameters. It is 
reasonable to expect the observable intensity to scale with the 
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Table 3. The values of the best fit parameters of the WS (sinu¬ 
soid with amplitude bo & phase (p) to the longitude variation 
(0° < I < 360°) of the statistical indicators SI of Galactic 
midplane, viz., the mean & median for : MSX PSC stars 
(|b| < 2.75°), intensity weighted & median latitude of emis¬ 
sion for MSX and DIRBE bands (|b| < 5°). The free parameter, 
k, for the null hypothesis is also tabulated. The ratio of per 
degree of freedom, R, for the best fitting null case and the WS 
case, is also presented (R> 1 indicates the sinusoid to be a bet¬ 
ter fit). 


A 

bo 

0 

k 

R 

jam 

(deg) 

(deg) 

(deg) 


MSX PSC : Mean, Median 

8.3 

0.18 ±0.02 

7.3 ±5.9 

-0.01 ± 0.02 

3.75 


0.28 ± 0.03 

13 ±7.3 

-0.02 ± 0.04 

2.78 

12.1 

0.22 ± 0.03 

1.3 ±8.5 

-0.05 ± 0.03 

2.15 


0.31 ±0.05 

1.1 ±9.2 

-0.07 ± 0.05 

1.98 

14.7 

0.24 ± 0.03 

1.9 ±7.8 

-0.04 ± 0.04 

2.46 


0.35 ± 0.06 

6.0 ± 9.7 

-0.05 ± 0.06 

1.95 

21.3 

0.28 ± 0.04 

18 ±9.0 

0.01 ± 0.04 

2.16 


0.43 ± 0.08 

32± 10 

0.07 ± 0.07 

1.83 

MSX Diffuse emission : Intensity weighted, Median 

8.3 

0.35 ± 0.05 

1.8 ±8.9 

0.02 ± 0.06 

2.19 


0.19 ±0.04 

-2.7 ± 13 

0.11 ±0.03 

1.18 

12.1 

0.43 ± 0.06 

-15 ±7.6 

0.01 ± 0.06 

2.64 


0.38 ± 0.06 

-23 ± 8.3 

0.06 ± 0.06 

2.30 

14.7 

0.58 ±0.10 

-12 ±9.6 

-0.03 ±0.10 

2.01 


0.99 ± 0.22 

-25 ± 13 

-0.16 ±0.19 

1.51 

21.3 

0.51 ±0.09 

-14± 10 

0.08 ± 0.08 

1.82 


0.54 ± 0.08 

-25 ± 8.7 

0.05 ± 0.08 

2.22 


DIRBE emission 

: Intensity weighted, Median 


1.25 

0.22 ± 0.07 

68± 17 

-0.08 ± 0.05 

1.18 


0.39 ±0.11 

65± 17 

-0.24 ± 0.08 

1.05 

2.2 

0.20 ± 0.05 

34± 14 

-0.04 ± 0.04 

1.41 


0.30 ± 0.08 

30± 15 

-0.14 ±0.06 

1.20 

3.5 

0.26 ± 0.04 

15 ±9.8 

-0.03 ± 0.04 

1.95 


0.37 ± 0.07 

11 ± 11 

-0.14 ±0.06 

1.55 

4.9 

0.30 ± 0.04 

6.8 ± 8.6 

-0.03 ± 0.05 

2.27 


0.41 ± 0.07 

3.8 ±9.3 

-0.14 ±0.06 

1.80 

12 

0.39 ± 0.05 

-13 ±7.0 

-0.04 ± 0.06 

2.89 


0.51 ±0.08 

-11 ±8.9 

-0.19 ±0.08 

1.85 

25 

0.35 ± 0.04 

-6.8 ±7.3 

-0.01 ± 0.05 

2.79 


0.49 ± 0.08 

-3.8 ±9.3 

-0.16 ±0.07 

1.86 

60 

0.54 ± 0.07 

-3.8 ±7.1 

-0.03 ± 0.08 

2.89 


0.68 ±0.11 

-6.1 ±9.6 

0.18 ±0.11 

1.87 

100 

0.49 ± 0.06 

-8.2 ±7.0 

-0.04 ± 0.07 

2.89 


0.64 ± 0.09 

-9.1 ±9.0 

-0.19 ±0.10 

1.95 

140 

0.52 ± 0.06 

-11 ±7.2 

-0.05 ± 0.08 

2.78 


0.65 ±0.10 

-11 ±8.9 

-0.20 ±0.10 

1.95 

240 

0.54 ± 0.07 

-13 ±7.4 

-0.04 ± 0.08 

2.72 


0.65 ±0.10 

-11 ±8.9 

-0.19 ±0.10 

1.97 


integral of the emission along the line of sight, p (Hammersley 
et al.[l995}. Hence, we write : 


rP~AUb) „ 

I{b) oc I e '!• '‘i dp (6) 

Jo 


for a given /, where R = -y + (p cos b)^ - 2 Rq p cos b cos I 
is the galactocentric distance and PmaxiUb) is the maximum 
length along a given line of sight which is only constrained 
by the finite size of the disk. These synthetic latitude profiles 
are found to be better represented by exponential functions than 
gaussian as expected. We have fitted exponential functions of 
the form poC *!< + rj io these latitude profiles, with po, iph and 
T] as free parameters. The tf/i, corresponds to characteristic an¬ 
gular scale and 77 represents the background (no latitude depen¬ 
dent background is anticipated due to the symmetry of the ge¬ 
ometry). The study has been carried out at two representative 
longitudes 30° and 90°. The lines of sight closer towards the 
Galactic centre have been avoided due to difficulty in compar¬ 
ing the calculations with measurements, since our description 
neglects the more complex details (like the bulge) which could 
be important. The following values of Ri : 1.5, 2.0, 2.5,..., 7.0, 
7.5 and 00 kpc (R/ = 00 corresponds to constant radial density 
distribution) have been explored. The set of Zh values explored 
are 40, 50, 60, ..., 90, 100, 125, 150,..., 275, 300 pc. 

The results of our calculations (for 1=30° and 90°) con¬ 
necting the characteristic angular scale (ye) with linear scale 
height (zh), for selected values of the radial scale lengths (R;), 
have been presented in Figure |3 The larger separation among 
curves corresponding to different values of R/ for /=90° as com¬ 
pared to 1=30°, is entirely expected from the shape and size 
of the Galactic disk and the heliocentric location of the ob¬ 
server. Overplotted as horizontal lines in this figure, are the 
values of y^ for mid infrared wavebands at corresponding lon¬ 
gitudes, as obtained earlier from fitting the observed latitude 
profiles (Section 3.2.2 and 3.3) for the 4 MSX bands and 2 
DIRBE bands (12 and 25 pm). This comparison helps con¬ 
straining the values of R/ and Zh for the distribution of warm in¬ 
terstellar dust. From the overlap of the observed and the model 
curves for I - 30°, we infer that the preferred range of Zh, is 
~ 60 - 100 pc, while R/ cannot be constrained due to the insen¬ 
sitivity to this parameter. However, the situation at I - 90° is 
quite favourable and provides a handle for R;. Assuming Zh to 
be in the above range, it is evident that the radial scale length, 
R;, should be less than ~ 6 kpc. Although the value of y^. ob¬ 
tained from DIRBE 25 pm band, 1.4°, seems incompatible with 
the above (Fig. 8 ,1 - 90°), we note that the values correspond¬ 
ing to the immediate neighbouring longitude bins are systemat¬ 
ically higher ( 2 . 0 ° & 2 . 2 °), consistent with the above constraint 
on R;. Giving equal weightage to the six mid-IR bands, we find 
the most favoured value of R; to be 2.81 + 0.63 kpc correspond¬ 
ing to Zh - 80 pc (central value of our implied range). Similarly, 
for the far infrared bands representing the cold interstellar dust, 
the scale heights Zh have been found to be ~ 60 pc. 

In addition, it is instructive to compare the characteristic 
angular scales, y^., obtained in the present study as a function 
of wavelength. For this, we have selected a moderate longitude 
bin (/ = 30°) where the complexities due to additional compo¬ 
nents of the Galactic structure near the Galactic center is ab¬ 
sent, yet the signal to noise ratio is reasonably good for reliable 
determination of y^s. The variation of je with wavelength, at 
1=30°, corresponding to the DIRBE and MSX wavebands (A > 
2.2 pm) have been presented in Figure |9l We find that y^. drops 
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with wavelength from near to mid infrared and remains more or 
less constant (within errors) thereafter upto far infrared wave¬ 
lengths. 

While the near infrared emission is expected to be domi¬ 
nated by stars, the mid infrared bands sample the warm inter¬ 
stellar dust in the vicinity of star forming regions. On the other 
hand, the far infrared emission is predominantly from cold dust 
heated by average interstellar radiation held. The observed de¬ 
crease in je obtained from 2.2 to 8.3 ;um bands could possibly 
be linked to the extinction due to dust. Typical wavelength de¬ 
pendence of dust extinction being oc /!““ (1 < a < 2), the ex¬ 
tinction is expected to be much larger for the near infrared as 
compared to mid or far infrared wavebands. The effect of ex¬ 
tinction is largest on the Galactic plane ib-Q) and progressively 
reduces at higher latitudes. The observed peak emission being 
reduced due to attenuation, the ht to the data could indicate an 
effectively larger angular scale height. 

5. Comparison with other resuits 

In this section we compare our results about the structure and 
WS of the Galactic disk with other studies of the distribution 
of stellar and the interstellar components. 

5.1. Indicators of structure of the Galactic disk 

The present study has constrained the scale height and ra¬ 
dial scale length for the distribution of warm interstellar dust 
(zh = 60-100 pc, Ri < 6 kpc; Section 4). Our Zh is close to that 
for the ISM (140 pc) from the successful model of the Galactic 
structure by Robin et al. (I2003> . although their R; is higher (4.5 
kpc). The distribution of dust studied by Freudenreich J1998t 
using near infrared measurements imply : Zh ~ 150 - 200 pc 
& Ri ~ 3.0 - 3.3 kpc, which is in reasonable agreement with 
our Ri. Drimmel & Spergel Il2001t . have obtained the values 
Z;,=134 pc & Ri-2.26 kpc for the cold dust component, using 
DIRBE 240 jum data, which compares very well with our re¬ 
sults. 

5.2. Signatures of Galactic warp 

In this subsection, our results of WS have been compared with 
other studies for the stellar, the warm and cold interstellar dust 
components. The Table 0]presents the summary of this overall 
comparison. 

5.2.1. Warp signature from the stellar component 

The WS for the stellar component obtained here using the 
MSX PSC, is in very good agreement with that obtained by 
Djorgovski & Sosin J1989t from the distribution of IRAS point 
sources, despite their sample covering a much larger latitude 
range of |b| < 10°. The estimates for the WS phase, cf), ob¬ 
tained here from both the approaches (using latitude profiles 
for the inner Galactic plane, LP; & using statistical indicators 
for the entire Galactic plane, SI) match extremely well with the 
IRAS results (see Table 0}. The WS amplitude obtained for 
the SI case is also consistent with the IRAS result, although the 


LP case (representing only the inner Galaxy) implies somewhat 
larger value. Similar consistency is also observed for the results 
obtained here using DIRBE data in L (3.5 pm) & M (4.9 pm) 
bands for both LP & SI. Lopez-Corredoira et al. J2002bt have 
analysed the data from Two Micron All Sky Survey (2MASS) 
in the context of Galactic warp. Their best fit warped model 
(Pig. 15, b - -1-3° lb- -3°) corresponds to a value of -i-2° 
for our phase parameter (p, which is again similar to those ob¬ 
tained here from all the four bands of MSX PSC as well as L & 
M bands of DIRBE. Preudenreich J1998> from a study of the 
Galactic disk obtain -0.08° < <p < 0.8°, which again is consis¬ 
tent with our results. 

It is interesting to note that, stellar distributions obtained 
from MSX and IRAS show fair amount of symmetry between 
the positive and negative peaks of the warp signature, unlike 
the case for H I distribution (Burton llOSSk 

5.2.2. Warp signature from the Interstellar dust 
component 

The diffuse mid infrared (8 -21//m) emission detected by MSX 
is expected to originate from the warmer (T ^ 100 K) compo¬ 
nent of the interstellar dust. In order to be heated to such high 
temperatures, the dust grains need strong local radiation fields, 
e.g. near young stellar objects (OB stars). Hence, detection of 
WS from diffuse mid infrared emission provides information 
about the distribution of warm dust around luminous young 
stars in the Galactic plane. Our results of the WS for the diffuse 
component obtained from the 4 bands of MSX and the two mid 
infrared bands of DIRBE (12 & 25 pm) have also been com¬ 
pared in Table 0] The cps obtained are consistent within errors. 
The amplitudes inferred from the MSX and DIRBE are similar 
with the former being slightly larger. An identical systematic 
shift (~ 20°) in <p can be noticed between the values obtained 
for the inner Galactic plane (from LP) and those for the entire 
disk (from SI), for all the mid infrared bands of MSX as well 
as DIRBE. 

The far infrared emission originates from the colder (T % 
30 K) component of interstellar dust in the Galactic plane. The 
WS detected in 60-240 pm wavebands of DIRBE should corre¬ 
spond to the distribution of this component. The 0s agree very 
well with those for the stellar as well as the warm dust com¬ 
ponent (Table 1^. The amplitude of the WS is similar to that 
obtained for the warmer dust. 

6. Summary 

Taking advantage of the recent high resolution (~ 18") Galactic 
plane survey in the mid infrared (8, 12, 14 & 21 pm) by the 
Midcourse Space Experiment (MSX), the distribution of stars 
(from the Point Source Catalog) as well as warm interstellar 
dust (from the panoramic images) in the Galactic disk have 
been investigated. The location of the Galactic mid-plane and 
the characteristic angular scales along latitude have been ob¬ 
tained from latitude profiles (-90° < I < 90°), and statistical 
indicators (entire Galactic plane). Using a simple geometric 
description of the Galactic disk, the angular scales have been 
translated to scale heights {zh) and radial length scales (R/). The 














8 


S. Vig et al.: Warp of the Galactic disk as seen in mid infrared 


Table 4. Comparison of parameters of warp signature (amplitude and phase) for the Galactic disk obtained from different studies. 


/l(//m) 

(Band) 

Survey 

Comp. 

LP/SF 

WS Amplitude 
(deg) 

WS Phase 
(deg) 

T2^ 

2MASS 

Stellar 



2 

(NIR) 






2.2° 

DIRBE 

Stellar 

LP 

0.5 

19 

(NIR) 



SI 

0.2 - 0.3 

30-34 

3.5, 4.9° 

DIRBE 

Stellar 

LP 

O 

1 

o 

15 

(NIR) 



SI 

0.3-0.4 

4-15 

8.3-21° 

MSX 

Stellar 

LP 

0.4-0.5 

3-7 

(MIR) 



SI 

0.2-0.4 

1-18 

12-60°' 

(MIR/EIR) 

IRAS 

Stellar 


0.2-0.3 

0-10 

00 

MSX 

Warm dust 

LP 

0.4 

18 

(MIR) 



SI 

0.2-0.4 

-3-2 

12-21° 

MSX 

Warm dust 

LP 

0.4-0.7 

-6-11 

(MIR) 



SI 

0.4-1 

-25--12 

12, 25° 

DIRBE 

Warm dust 

LP 

0.4 

11 - 18 

(MIR) 



SI 

0.4-0.5 

-13--4 

60-240° 

DIRBE 

Cold dust 

LP 

d 

1 

d 

12-16 

(FIR) 



SI 

0.5-0.7 

-13--4 


“ From Latitude Profiles, -90° < / <90° (LP) / Statistical Indicators, 0° < Z < 360°(SI) 
Lopez-Corredoira et al. <2002b> 

Present work 

Djorgovski & Sosin 11989> 


following constraints have been obtained for the distribution of 
warm (T ^ 100 K) interstellar dust in the Galactic disk : Rj < 6 
kpc and 60 ^ Zh 5 100 pc. These are consistent with results 
from other studies. 

Signatures of warp-like feature in the distribution of stars 
as well as warm interstellar dust have been searched for using 
the inferred locations of the Galactic mid-plane as a function 
of longitude. The warp signature (WS) has been detected in all 
the cases : for point sources as well as diffuse emission in each 
of the 4 MSX bands. The WS for each case has been quantified 
by two parameters depicting its amplitude and phase (longitude 
of the line of nodes). 

Our analysis/study of the mid infrared emission has been 
extended to cover the entire infrared waveband by using the 
data from the DIRBE/COBE experiment (10 bands; 1.25-240 
Despite limited angular resolution of DIRBE, signatures 
of warp have been detected in each band with very similar 
phase as found from the MSX data. The inferred amplitudes 
of WS lie in the range : 0.2° - 0.7° (for A > 2.2//m) in lati¬ 
tude. Results from the present study have also been compared 
with those from other similar studies based on the IRAS and 
2 MASS surveys. 
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Ldtitude (degrees) 
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Fig. 1. The number density of stars as a function of latitude (lat¬ 
itude profile) detected in the A band of MSX PSC for the lon¬ 
gitude bins 30°, 60° and 180°. The solid and the dotted lines 
represent best fit exponential and gaussian functions, respec¬ 
tively (for 1= 30° & 60°). 



0 100 200 300 
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0 100 200 300 

Longitude (degrees) 




Fig. 2. The variation of the location of Galactic midplane with 
longitude, /1(0, as obtained from fitting exponential to the lat¬ 
itude profiles (LP) for our MSX PSC sample (histogram) for 
the inner Galaxy in 4 MSX bands. The solid lines refer to the 
extracted WS in (a) 8 pm, (b) 12 pm, (c) 14 pm, and (d) 21 pm 
bands. 
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Longitude (degrees) Longitude (degrees) 

Fig. 3. The variation of the Sis of the Galactic midplane, viz. 
mean (thick line histogram) and median (thin line histogram) 
latitude as a function of longitude for our MSX PSC sample in 
the 4 MSX bands : (a) 8 jjm, (b) 12 jsm, (c) 14 jim, and (d) 21 
yum. The solid and the dotted curves represent the extracted WS 
from the mean and median latitudes, respectively. 
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Fig. 5. The variation of the location of Galactic midplane with 
longitude as obtained by fitting exponential to the latitude pro¬ 
file (LP) of the total (diffuse) emission measured in 4 MSX 
bands (histogram) for the inner Galaxy. The solid lines refer to 
the extracted WS in (a) 8 jim, (b) 12 jim, (c) 14 //m, and (d) 21 
jum bands. 



Lotitude (degrees) 


Fig. 4. Latitude profiles of MSX A band intensity for longi¬ 
tude bins 10°, 60° and 180°. The solid and the dotted lines rep¬ 
resent best fitting exponential and gaussian functions, respec¬ 
tively (for f=10° & 60°). 
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Longitude (degrees) Longitude (degrees) 

Fig. 6. The variation of the Sis of the Galactic midplane, viz. in¬ 
tensity weighted mean {Bintwt', thick line histogram) and median 
{Bmed', thin line histogram) latitudes as a function of longitude 
for the diffuse emission in the 4 MSX bands : (a) 8 yum, (b) 12 
yum, (c) 14 yum, and (d) 21 yum. The solid and the dotted curves 
refer to the extracted WSs for the intensity weighted and me¬ 
dian latitudes, respectively. 
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Fig. 8 . Relations between the linear (z/,) and the angular (y^.) 
scale heights from a geometric description of the Galactic disk 
(for a heliocentric observer), for different assumed radial scale 
lengths, R;. The two panels represent longitudes I = 30° & I - 
90°. The y^s obtained from the exponential fits for the diffuse 
mid infrared emission in 4 MSX bands (for A, C, D & E bands 
: 0.85°, 0.92°, 0.96° & 0.67° for /=30°; and 3.08°, 2.45°, 2.43° 
& 2.37° for /=90°) and the two DIRBE bands (for 12 & 25 
bands ; 0.84°, 0.72° for 1=30°; and 2.81°, 1.40° for /=90°) are 
also plotted as horizontal solid lines for comparison. 


Fig. 7. Same as in Fig ISlbut for DIRBE intensity data in the ten 
bands: (a) 1.25 //m,(b) 2.2 fim, (c) 3.5 jum, (d) 4.9 /zm, (e) 12 
//m, (f) 25 jum, (g) 60 /rm, (h) 100 jum, (i) 140 /rm and (j) 240 
//m band. The solid and the dotted lines refer to the extracted 
warp signatures corresponding to the best fit sinusoid to 
and Bmed^i respectively. 
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Fig. 9. Wavelength dependence of the angular scale height, 
for the longitude bin 1=30°. The open circles represent the 
DIRBE wavebands for A > 2.2 /rm and the solid squares repre¬ 
sent the diffuse emission in the MSX bands at 8, 12, 14 and 21 
yum. The error bars indicate Icr uncertainity. At longer wave¬ 
lengths the errors are smaller than the symbol size. 



